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Abstract 

The charged multiplicity distributions and the mean charged multiplicity have 
been investigated in inclusive neutral current deep inelastic ep scattering with 
the ZEUS detector at HERA, using an integrated luminosity of 38.6 pb" 1 . The 
measurements were performed in the current region of the Breit frame, as well as 
in the current fragmentation region of the hadronic centre-of-mass frame. The 
KNO-scaling properties of the data were investigated and the energy dependence 
was studied using different energy scales. The data are compared to results 
obtained in e + e~ collisions and to previous DIS measurements as well as to 
leading-logarithm parton-shower Monte Carlo predictions. 
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1 Introduction 



The production of multi-hadronic final states in high-energy two-body collisions has long 
been a subject of great interest from experimental and theoretical points of view. The 
charged multiplicity at HERA has been measured previously by the HI [1-4] and ZEUS [5- 
7] experiments. In this paper, new measurements by the ZEUS collaboration of the 
charged multiplicity in deep inelastic scattering (DIS) are presented. 

Measurements are performed in the hadronic centre-of-mass (HCM) frame and the results 
are compared with those obtained in e + e~ collisions, as well as with those from previous 
DIS experiments [1,8-10]. For the ep final state, differences are expected in the photon 
(current) and proton (target) fragmentation regions, due to the asymmetric nature of the 
reaction. The detector acceptance only allows the current fragmentation region to be 
measured. 

Measurements of the charged multiplicity are also performed in the current region of the 
Breit frame, which should behave similarly to one hemisphere in e + e~ collisions. Previous 
DIS results [2,5, 11] using Q, the virtuality of the exchanged photon, as the scale showed 
a reasonable agreement with e + e~ data. However, this agreement degraded at values of 
Q below 6-8 GeV. In this paper, the energy of the current region of the Breit frame is 
used as the scale to compare with e + e~ data. 

An alternative energy scale, the invariant mass of the hadronic system, has also been used 
in both the Breit and HCM frames. The results using this variable are also compared to 
results from e + e~ collisions. 



2 Experimental set-up 

The data were collected with the ZEUS detector during the 1996 and 1997 running periods, 
when HERA operated with protons of energy E p = 820 GeV and positrons of energy 
E e = 27.5 GeV, and correspond to an integrated luminosity of 38.6 ± 0.6 pb" 1 . 

The ZEUS detector is described in detail elsewhere [12]. The most important components 
used in the current analysis were the central tracking detector and the calorimeter. 

Charged particles were tracked in the central tracking detector (CTD) [13], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consisted 
of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the polar-angle§ 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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region 15° < 6 < 164°. The transverse-momentum resolution for full-length tracks was 
°{Pt)/pt = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [14] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a(E)/E = O.lS/y/E for electrons and 
a{E)/E = 0.35/ VE for hadrons, with E in GeV. 

3 Data selection 

Deep inelastic scattering events were selected by requiring that the outgoing positron was 
measured in the CAL. The scattered-positron identification was based on a neural-network 
algorithm using the CAL information [15]. 

For the reconstruction of the photon virtuality, Q 2 , Bjorken x, and the 7*P centre-of-mass 
energy, W, the double angle (DA) method was chosen, in which the scattered-positron 
angle, 9 e , and the hadronic angle 7# are used [16]. In the naive quark-parton model, 7# 
is the angle of the scattered massless quark in the laboratory frame. 

For each event, the measurement of the charged multiplicity was performed using tracks 
reconstructed in the CTD. The energy of the hadronic final state was measured using a 
combination of track and CAL information, excluding the cells and the track associated 
with the scattered positron. The selected tracks and CAL clusters were treated as massless 
Energy Flow Objects (EFOs) [17]. The clustering of objects was done according to the 
Snowmass convention [18]. The transverse momentum, p?, of each EFO was required to 
be greater than 0.15 GeV. 

The event selection criteria were: 

• E e i > 12 GeV, where E e < is the energy of the scattered positron, to select neutral 
current DIS events; 

• He < 0.95, where y e is the scaling variable y as determined from the energy and polar 
angle of the scattered positron. This cut reduces the photoproduction background; 

• Z/jb > 0.04, where i/jb is the estimate of ?/ using the Jacquet-Blondel (JB) method [19]. 
It is defined as yjB = Ylihi^h — Pz h )/2E e , where the sum runs over all EFOs and E h 
and Pz h are the energies and longitudinal momenta of the EFOs respectively. This 
requirement guarantees sufficient accuracy for the DA reconstruction method; 



2 



• 35<<5<60GeV, where S = ~ -PzJ an d the sum runs over all EFOs and the 
scattered positron, to remove photoproduction events and events with large radiative 
corrections; 

• R > 25 cm, where R is the distance from the beam axis to the impact position of the 
scattered positron on the face of the CAL. This ensured that the positron was fully 
contained within the detector and its position reconstructed with sufficient accuracy; 

• \Z vtx \ < 50 cm, where Z vtx is the longitudinal position of the vertex, to reduce back- 
ground events from non-ep collisions. 

To ensure high-quality tracks reconstructed with high efficiency, the following require- 
ments were made: 

• the tracks had to pass through at least three CTD superlayers; 

• the tracks had to be associated with the primary event vertex; 

• the tracks were restricted to the region |?7lab| < 1.75, where 7/lab is the pseudorapidity 
of the tracks in the laboratory frame; 

• the tracks had pr > 0.15 GeV. 

The analysis was restricted to the kinematic range Q 2 > 25 GeV 2 and 70 < W < 225 GeV. 

4 Analysis method 

Due to the large asymmetry of the beam energies at HERA, a large fraction of the hadronic 
final state close to the proton direction lies outside the detector acceptance. Therefore 
only hadrons belonging to the current fragmentation regions of the HCM and Breit frames 
were used in this analysis. The boost to the corresponding reference frames was calculated 
using the positron four-momentum taken from the DA method. 

4.1 Breit frame 

In the Breit frame, which is defined by the condition that the momentum of the exchanged 
virtual boson is purely spacelike, q = (0, 0, —Q), the particles produced in the interaction 
can be assigned to one of two regions: the current region if their longitudinal momentum 
in the Breit frame is negative, and the target region if their longitudinal momentum is 
positive. The hadronic system of the current region used in this analysis is almost fully 
(about 95%) contained within the acceptance of the CTD. 
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Previous analyses compared the mean charged multiplicity (n^) as a function of Q to 
(n c h)/2 as a function of \fs in e + e~ collisions [2,5]. For values of Q > 8 GeV, reason- 
able agreement was observed, while some disagreement was found for Q < 8 GeV. The 
difference can be understood in terms of higher-order processes [20], which change the 
available energy in the current region of the Breit frame, E^, which is no longer equal to 
Q/2. In this analysis the quantity 2 ■ E# is used as a scale. On an event-by-event basis 
this method should compensate for particles and their corresponding energies migrating 
between current and target regions of the Breit frame. 

4.2 Hadronic centre-of-mass frame 

In the HCM frame, the exchanged virtual boson has four-momentum q = {E,p} = 
( W 2 w* j 0; 0; ~~2W )• hadronic final state is separated into the photon (current) and 
proton (target) fragmentation regions. About 60 - 80% of the current region of the HCM 
frame is contained within the acceptance of the CTD. 

The multiplicity in the HCM frame in DIS is usually studied as a function of W [1,8-10]. 
At HERA, the energy in the current region of the HCM frame, -&hcm> coincides with H^/2 
within 0.3 — 0.4%. Thus, for practical reasons, W was used as the energy scale. 

4.3 Invariant mass of the hadronic system 

Charged multiplicities in the current region of both the Breit and HCM frames were also 
measured as functions of the invariant mass: 

where the sum runs over all charged and neutral particles of the corresponding hadronic 
system. 

5 Monte Carlo models, acceptance corrections and 
systematic errors 

Samples of neutral current DIS events were generated using the colour-dipole model as 
implemented in Ariadne 4.12 [21] or with the MEPS model of Lepto 6.5 [22]. Both 
programs were interfaced to Heracles 4.6.1 [23] using the Djangoh 1.1 [24] program. 
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Both Ariadne and Lepto use the Lund string model [25] for the hadronisation. The 
event samples were generated using the CTEQ4D [26] parameterisation of the parton 
distribution function (PDF) in the proton. The cluster hadronisation model as imple- 
mented in Herwig 6.1 [27] was used to estimate the effect of different hadronisation 
schemes on the unfolding procedure. The minimum transverse momentum of outgoing 
partons in the hard interaction and of partons participating in multi-parton interactions 
is adjustable in Herwig using the parameter p™ m . The parameter was tuned in this 
analysis to improve the description of the detector distributions. The best agreement was 
found for pf n = 2.5 GeV [28]. The HERWIG samples were generated using CTEQ5L [29] 
proton PDF parametrisation. The MC samples were used both for data correction and 
for comparison of the data to the model predictions. 

The corrections applied to the data accounted for the effects of acceptance and resolution 
of the detector, event selection cuts, QED-radiative effects, track reconstruction, track 
selection cuts, and energy losses in the inactive material in front of the calorimeter in 
the case of the energy measurement. Finally the multiplicity distributions were corrected 
using a matrix unfolding method as described in earlier studies [5] . 

The generated events were passed through a full simulation of the detector, using a pro- 
gram based on Geant 3.13 [30], and processed and selected with the same programs as 
used for the data. The simulated samples were used to determine the response of the 
detector and to evaluate the correction factors necessary to obtain hadron-level quanti- 
ties. The hadron level is defined by those hadrons with lifetime r > 3 • 10 -10 seconds. 
In order to compare the results to different experiments, corrections were calculated both 
including and not including the decay products of Kg and A. 

The dominant systematic uncertainty in this analysis arises from the simulation of the 
hadronic final state. To correct the data, the average of the correction factors from the 
Ariadne and Herwig MC programs was used. One half of the difference, as large as 
5%, was assigned to the systematic uncertainties. 

Other sources of uncertainty are (typical values of the uncertainties are shown in paren- 
theses): event reconstruction and selection (< 0.5%), track reconstruction and selection 
(< 0.5%), and the uncertainty due to variation of the Q 2 cut by its resolution (< 1.7%). 
The uncertainty due to contamination from diffractive events is negligible. The individ- 
ual systematic uncertainties were added in quadrature. The major correlated uncertainty 
comes from the CAL energy scale (< 1.5%) and is not shown in tables and figures. A 
detailed study of the sources contributing to the uncertainties of the measurement can be 
found in [31]. 
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6 Results 



6.1 Multiplicity distributions 

The multiplicity distributions in the current region of the Breit frame are presented in 
Fig. [1] and Table [1] and in Fig. [2] and Table [2] in bins of 2 • E# and in bins of M e g 
respectively. The kinematic range of the analysis restricts the M e g measurements in the 
current region of the Breit frame to a maximum value of about 20 GeV. The predictions of 
Ariadne, Lepto and Herwig are also shown. All three MC models generally describe 
the data, but Ariadne gives the best description. In all Figures and Tables the charged- 
particle decay products of Kg and A are included, unless otherwise stated. 

For a given bin of energy, the multiplicities as a function of 2 • E# and M c g differ by 
approximately a factor of two. This is due to the fact that 2 ■ (as well as Q or W) 
characterise the total centre-of-mass energy of the system of which only one hemisphere 
is measured. On the other hand the M e g method measures the multiplicity of the system 
with respect to the corresponding invariant mass. 

The multiplicity distributions in the current region of the HCM frame are presented 
in Fig. [3] and Table [3] in bins of W. Both Ariadne and Lepto predict similar W 
distributions and give a reasonable description of the data. Herwig predicts longer tails 
for the multiplicity distributions in bins of W. This leads to higher average multiplicities, 
affecting the unfolded multiplicity values and increasing the systematic uncertainties of 
the measurement. The Ariadne predictions vary slightly from both Lepto and Herwig 
at the rising edge of the distributions, which also leads to an increase of the systematic 
uncertainties. 

The multiplicity distributions in the current region of the HCM frame in bins of M e g are 
presented in Fig. H] and Table HI To minimise the extrapolation both in multiplicity and 
Mgg, an additional requirement, |?7lab| < 1-75, was applied at the hadron level. None of 
the MC models shown in Fig. H] give a complete description of the data. This is most 
visible at higher values of M e g. 

6.2 KNO scaling 

The multiplicity distributions are expected to scale with energy as discussed in detail 
elsewhere [32]. A commonly used form of the scaling, from KNO [33], is shown in Figs. 5- 
8, where the product of the multiplicity distribution P{n c h) with average multiplicity 
(n c h), (^ch)-P(^ch), is shown as a function of n c h/{n c h). 

In Fig. 0(a), the KNO distributions measured in bins of W in the current region of the 
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HCM frame are shown. Within the uncertainties, the distributions measured in three 
bins of W agree. They also agree well with the average distribution, which was calculated 
using data for the entire W region, 70 < W < 225 GeV. This average KNO spectrum, 
presented as a histogram, is shown in Figs. [5] and [6] as a reference KNO distribution. 

The reference KNO distribution is compared to the measurements in the current region of 
the Breit frame in Figs. [3(b) and [5(c) in bins of 2 ■ For values of 2 • E% > 12 GeV, as 
shown in Fig. 0(c), the measurements are in reasonable agreement with the reference KNO 
histogram. For lower values of 2 • E^, as demonstrated in Fig. [5](b), the distributions do 
not follow the KNO-scaling behaviour; they have different shapes, but approach the KNO 
curve with increasing values of 2 • E#. 

In Fig.[6j the KNO distributions measured in bins of M e g in the current region of the HCM 
frame are presented. The multiplicity distributions in bins of M c $ do not follow the same 
KNO scaling observed for measurements as functions of W or 2 ■ but do demonstrate 
scaling behaviour for M e s values above 8 GeV. The measurements at M e s < 4 GeV 
in the current regions of both the Breit and HCM frames behave differently from the 
measurements at the higher values of M e g. 

Figure [7J shows a comparison of the KNO distributions from ZEUS with results obtained 
in e + e~ collisions. The measurements in bins of 2 • E^, for 2 • E# > 12 GeV, and in 
bins of W are plotted together and compared with measurements in one hemisphere of 
e + e~. In Fig. [7(a), a comparison with results from the TASSO collaboration [34] in the 
energy range 14 < Js^ e < 44 GeV is shown. At LEP only DELPHI [35] and OPAL [36] 
performed measurements in a single hemisphere at J~s^e = 91.2 GeV. A comparison with 
the present data is shown in Fig. [7Jb). The systematic uncertainties are not shown, but 
within statistical uncertainties there is a remarkable agreement between ep and e + e~ 
results. However, the LEP data differ somewhat from the present measurement in the 
peak region and at very low values of n c h/ (^ c h)- 

The data as a function of M e g, for M e g- > 8 GeV, are compared with the e + e~ measure- 
ments for the whole event in Fig. El Both the TASSO and LEP [35-37] data (91.2 < 
a/s^ < 209 GeV) agree with the present measurement. 

6.3 Mean charged multiplicity 

Figure [9] and Table [5] show the mean charged multiplicity, (n c h), in the current region 
of the HCM frame as a function of W and the mean charged multiplicity in the current 
region of the Breit frame as a function of 2 • E%. The Kg and A hadrons were considered 
stable in Fig. [91 where the data are compared with results of previously published HERA 
measurements [1,2,4,5]. As expected, at low values of 2 ■ E#, the measurement differs 
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with those as a function of Q (see Section l4TTj) . At higher values of 2 • Eg the data agree 
within the experimental uncertainties with the previous ZEUS and HI measurements, 
but lie systematically above them. The data are in good agreement with the Ariadne 
and Lepto predictions. The Herwig predictions also describe the data but are below 
those from Ariadne and Lepto. In the current region of the HCM, the measurement 
agrees well, with improved statistical and systematic uncertainties, with the earlier HI 
results. The Ariadne and Lepto predictions agree with the data. Herwig predicts a 
very different slope, with much higher multiplicities at higher energies; with increasing 
energy the agreement with data degrades. 

The mean charged multiplicities in the current regions of the Breit and HCM frames are 
presented in Fig. [10] and Table [6] as a function of the invariant mass of the corresponding 
hadronic system, M eff . In Figs. [IDT a) and [TOTb). the multiplicities are compared to the 
MC predictions. All three MC models describe the data reasonably well; however in 
the last M e fj bin in the current region of the HCM, the Herwig prediction is too high. 
In Fig. [TOT c). both measurements are shown together and compared with different MC 
curves calculated using the Ariadne MC. The measurements in the Breit and HCM 
frames agree at values of M e g less than 10 GeV. Above this value, (n c h) rises much faster 
with M e ff in the current region of the HCM frame than in the current region of the Breit 
frame. Since the HCM measurement was restricted in 77, a separate Ariadne calculation 
was performed in the total current region of the HCM frame. The difference is small, 
although the rise of (n c h) with M e s is faster in the total current region of the HCM frame. 

Figure [TOTc) also shows 2 • (n c h) as a function of 2 • Eg. This measurement exhibits the 
same behaviour as (n c h) as a function of M e g in the current region of the Breit frame but 
differs from that in the HCM frame. The multiplicity in the current region of the HCM 
frame rises much faster with the invariant mass than with W. 

Finally, Fig. [TT] combines the mean charged multiplicities measured in the current regions 
of the Breit and HCM frames as functions of the respective energy scales, 2 ■ Eg and W. 
Also shown are the measurements from e + e~ [34-38] and fixed-target [8-10] experiments. 
The fixed-target data were scaled by a factor two, since they only measure one hemisphere 
and by a factor 1.08, to correct for the decays of the K s and A, as estimated using 
Ariadne MC. 

The measurements presented in this paper show good overall agreement with those from 
other experiments, exhibiting approximately the same dependence of the mean charged 
multiplicity on the respective energy scale. At low values of the energy, the measurement 
as a function of 2 ■ Eg agrees well with e + e~ data, in contrast to the previous measure- 
ments as a function of Q [2,5]. The measurements in the current region of the HCM 
agree with the LEP data, but are systematically below them. The data from fixed-target 
DIS experiments [1,8-10] deviate from the observed energy dependence at energies above 
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15 GeV. The Ariadne MC prediction generally describes the energy dependence of the 
data over the entire region. However, the prediction in the HCM frame is generally lower 
than the data and than the prediction in the Breit frame. The Herwig MC model does 
not give a good overall description of the data. 

7 Summary and conclusions 

The charged multiplicity distributions and the mean charged multiplicity have been in- 
vestigated in inclusive neutral current deep inelastic ep scattering in the kinematic range 
Q 2 > 25 GeV 2 and 70 < W < 225 GeV in terms of different energy scales. The scale 
2 ■ E B T , was used in the current region of the Breit frame. In the current region of the 
HCM frame, W was used and the invariant mass, M e g, was used in both frames. 

In terms of KNO scaling, the charged multiplicities in the current regions of the Breit and 
HCM frames exhibit the same behaviour as those in one hemisphere of e + e~ collisions 
when 2 • E B or W are considered. When the energy scale M c g is used, the charged 
multiplicities exhibit the same KNO-scaling behaviour as those for the whole e + e~ event. 

The mean charged multiplicity in the current region of the Breit frame scales with M e g 
in the same way as 2 ■ (n^) scales with 2 ■ E B and, therefore, as (n c h) scales with ^/s^ in 
e + e~ collisions. The mean charged multiplicity in the current region of HCM frame as a 
function of M e g rises faster than that in the current region of the Breit frame. 

The energy scale 2 • E B , rather than Q, gives better agreement between the mean charged 
multiplicity measured in the current region of the Breit frame and that measured in e + e~ 
collisions. The measurements of (n c h) as a function of W agree, within the uncertain- 
ties, with the data from e + e~ collisions. When using these scales, ep DIS data can be 
consistently compared with data from e + e~, fiP and vP scattering over a wide energy 
region. 
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Table 1: Multiplicity distributions 100 • l/NdN/dn c h measured in the current 
region of the Breit frame in bins of 2 ■ E# . The first errors are statistical and the 
second are the systematic uncertainties. 
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Table 2: Multiplicity distributions 100 • \/NdN/dn c h measured in the current 
region of the Breit frame in bins of M e ^. The first errors are statistical and the 
second are the systematic uncertainties. 



14 



W (GeV) 


70 - 100 


100 - 150 


150 - 225 


n c h=0 


0.09 ± 0.03 + nn} 

— U . U 1 


0.11 ± 0.03+?,' no 

— U.Uz 


0.09 ± 0.02 + nnt 

— U.U4 


1 


0.37 ± 0.04+!! )'i 

— U . 1 o 


0.26 ± 0.03+?,' )9 

— U. 1U 


0.25±0.04 + nl7 


2 


1.50 ±0.07+2 ™ 


1.02 ± 0.06+?, 1% 

— 0.2o 


0.81±0.05 + nl„ 
— u. iy 


3 


2.58 ± 0.09+noi 

— U.OO 


1.55 ± 0.06+n,7 

— U.l / 


0.99±0.05 + nl7 

— U.l/ 


4 


5.22 ± 0.13+n'S 


3.60 ± 0.09+n tl 

— u. by 


2.43±0.08 + nKfi 

— U . bo 


5 


6.47 ±0.14+S'?2 

— U.Jo 


4.43 ± 0.10+n J! 

— U . 1 o 


3.00 ± 0.08+2'?! 

— u.ub 


6 


8.97±0.16 + n'5o 

— U.oo 


6.71 ± 0.12+°'^ 

— U.85 


4.72±0.10 + non 
— 0. oU 


7 


9.57±0.17+?w 1 ,q 
— u.uy 


7.40±0.13+nls 

— U ■ J-O 


5.51 ± 

— U.^O 


8 


9.89 ± 0.16+n'II 

— U. (o 


8.41 ±0.13+°'^ 
— u.oy 


6.82 ± 0.12+2-SI 

— U.oY 


9 


9.53 ± 0.16+n? 3 , 

— U.oo 


8.45 ± 0.13+n'r} 

— U.Ol 


7.18 ± 0.12+2 K 

— U.o i 


10 


8.73 ± 0.16+n 19 

— U. (0 


8.40 ±0.13+?,'^ 

— U.oo 


7.61 ± 0.13+nIl 

— U. ( 4 


11 


7.82 ± 0.15±g;J} 


8.05±0.13lg;l? 


7.54±0.13 + n5^ 

— U.O i 


12 


6.63±0.14+! ) ;}f; 


7.31 ± 0.12+°-^ 
— u.4y 


7.41 ± 0.13+°-^ 

— U.oo 


13 


5.50 ± 0.13+nli 

— U. 12 


6.46 ± 0.12+nlS 

— U. lo 


6.82±0.12 + n 3 n 

— U.oU 


14 


4.32 ± 0.12+2-J? 

— U. 15 


5.53 ± 0.11+?' n? 

— 0.02 


6.33 ± 0.12 + ° |1 

— . 6 1 


15 


3.45 ± 0.11 + n?f 

— U.ZO 


4.72 ± 0.10+?.' lo 
— u.iy 


5.69 ±0.11+° ^ 

— U.Uo 


16 


2.58 ± 0.09+fJol 

— U.o i 


3.90 ±0.10+? :^ 

— 0.2 I 


4.93±0.11+0-?» 

— U.Uo 


17 


1.96 ± 0.09+n 11 

— U . 44 


3.16 ± 0.09+0'H 

— U.44 


4.22±0.10 + n3n 

— U . oU 


18 


1.46 ± 0.08+S-1S 

— U-45 


2.58 ± 0.08+?' 1? 

— U.4o 


3.48 ± 0.09 + n on 

— U . oU 


19 


1.05 ± 0.07+n« 

— U.oo 


1.98 ± 0.08±2JJ 


2.95±0.09 + no7 

— U.o l 


20 


0.74 ± 0.06+n?,l 

— U.o4 


1.50 ± 0.07+nil 

— U.Ol 


2.32 ± 0.08 + ?i5 

— U.4o 


21 


0.51 ± 0.06+n?K 

— U./o 


1.14±0.06 + nls 

— U.4o 


1.91±0.08 + n5? 

— U . o 2 


22 


0.34 ± 0.05+n'lq 


0.90±0.06 + n?s 

— U.oo 


1.51 ± 0.07+nll 

— U.4o 


23 


0.25 ± 0.05±g;];{j 


0.65 ± 0.05+?' It 


1.22 ±0.07+2-1; 

— U.4 ( 


24 


0.17 ±0.041^2 


0.46 ± 0.05±°^ 


0.94 ± 0.06i{!;ll 


25 


0.10±0.04+.°;°f 


0.36 ± 0.05lo;23 


0.74±0.06+.g;lg 


26 


0.08 ± 0.03lo Qg 


0.28 ± 0.04j;!5;2o 


0.58 ± 0.05lo Jl 


27 


0.05 ± 0.03lg'°l 


0.19 ± 0.041°;^ 


0.45 ± 0.05i°;27 


28 


0.03 ± 0.02i°;o3 


0.15 ± 0.04l°^2 


0.36 ± 0.05lo 23 


29 


0.02 ± 0.02^^2 


0.10 ± 0.04+.S/2J 


0.29 ± 0.05lo 20 


30 


0.01 ± 0.03to;° 2 


0.07 ± 0.02±g 


0.22±0.05i°;l7 


31 




0.05 ± 0.03+^1 


0.17±0.04i55^2 


32 




0.04 ± 0.03lg;°3 


0.13±0.05l°il 


33 




0.03 ± 0.03+.g;°3 


0.10 ±0.04+.°;°* 


34 




0.01 ± 0.02±g ;^ 


0.08 ± 0.03lo gg 


35 
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Table 3: Multiplicity distributions 100 • l/NdN/dn c h measured in the current 
region of the HCM in bins ofW. The first errors are statistical and the second are 
the systematic uncertainties. 
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Table 4: Multiplicity distributions 100 ■ l/NdN/dn c h measured in the current 
region of the HCM frame in bins of M c g . The first errors are statistical and the 
second are the systematic uncertainties. 
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2 ■ Eg (GeV) 



2.9 

5.9 

9.6 

14.8 

23.8 
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W (GeV) 
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68 ± 





03 
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Table 5: Mean charged multiplicity, (n^), measured in the current region of the 
Breit frame as a function of 2 • Eg and in the current fragmentation region of the 
H CM frame as a function ofW. The first errors are statistical and the second are 
the systematic uncertainties. 



Current region of the Breit frame 
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(71^) (if , A decay) 
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3.38 ± 0.02+°;°^ 
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10.84 ± 0.02+8;^ 
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14.17 ± 0.05±g;lg 



Table 6: Mean charged multiplicity, (n c h), measured in the current region of 
the Breit frame and in the current fragmentation region of the HCM frame as a 
function of M e g. The first errors are statistical and the second are the systematic 
uncertainties. 
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Figure 1: Multiplicity distributions measured in the current region of the Breit 
frame in bins of 2 ■ E§ (solid circles). The inner error bars represent the statistical 
uncertainties and the outer error bars the statistical and systematic uncertainties 
added in quadrature. The predictions (solid lines) of different MC models are also 
shown. 
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Figure 2: Multiplicity distributions measured in the current region of the Breit 
frame in bins of M e g (solid circles). The inner error bars represent the statistical 
uncertainties and the outer error bars the statistical and systematic uncertainties 
added in quadrature. The predictions (solid lines) of different MC models are also 
shown. 
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Figure 3: Multiplicity distributions measured in the current region of the HCM 
frame in bins of W (solid circles). The inner error bars represent the statistical 
uncertainties and the outer error bars the statistical and systematic uncertainties 
added in quadrature. The predictions (solid lines) of different MC models are also 
shown. 
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Figure 4: Multiplicity distributions measured in the current region of the HCM 
frame in bins of M e g (solid circles). The inner error bars represent the statistical 
uncertainties and the outer error bars the statistical and systematic uncertainties 
added in quadrature. The predictions (solid lines) of different MC models are also 
shown. 
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Figure 5: The multiplicity distributions plotted in KNO form. The reference KNO 
histogram represents the measured distribution in the HCM frame in the entire W 
range, (a) Multiplicity distributions measured in the current region of the HCM 
frame in bins ofW. Multiplicity distributions measured in the current region of the 
Breit frame in bins of 2 ■ for (b)1.5<2- < 12 GeV and (c) 12 < 2 • E% < 
100 GeV. 
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Figure 6: The multiplicity distributions plotted in KNO form and compared 
to the reference KNO distribution (histogram). The multiplicity distributions are 
measured in bins of M eS in the current regions of (a) the HCM frame, restricted in 
Vlab, (b) the Breit frame. 
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Figure 7: The multiplicity distributions plotted in KNO form and compared to the 
results of e + e~ experiments. The solid circles represent the ZEUS data measured 
in the current region of the Breit frame in bins of 2 • E# , for 2 • E# > 12 GeV, 
and the solid squares represent the data measured in the current region of the HCM 
frame in bins ofW. Multiplicities measured in one hemisphere of the e + e~ collision 
are shown in bins of \fs^ e for (a) the TASSO experiment [34], in the energy range 
14 < ^fs^e < 44 GeV, and (b) for the LEP experiments [35, 36], measured at energy 
^fsZ= 91.2 GeV. 
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Figure 8: The multiplicity distributions plotted in KNO form and compared to 
the results of the e + e~ experiments. The solid circles represent the ZEUS data 
measured in the current region of the Breit frame and the solid squares represent 
the data measured in the current region of the HCM frame both in bins of M e g, 
for M c g > 8 GeV. Charged multiplicities measured for the whole event in e + e~ 
collisions are shown in bins of ^fs^ e for (a) the TASSO experiment [34], in the 
energy range 14 < yfs^ e < 44 GeV, and (b) for the LEP experiments [35-37] in 
the energy range 91.2 < J~s[T e < 209 GeV. 
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Figure 9: Mean charged multiplicity, (n c h), in the current region of the Breit 
frame as a function of 2 • E^ and in the current fragmentation region of the HCM 
frame as a function ofW. The inner error bars represent the statistical uncertain- 
ties, typically smaller than the size of the symbols. The outer error bars represent 
the quadratic sum of statistical and systematic uncertainties. Also shown are the 
results of previous HERA measurements [1, 2,4, 5] and predictions from Ariadne, 
Lepto and Herwig. The decay products of Kg and A are not included. 
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Figure 10: Mean charged multiplicity, (n^) , measured as a function of M e g 
(a) in the current region of the Breit frame and (b) in the current region of the 
HCM frame compared to MC predictions, (c) Comparison between measurements 
in the current regions of the Breit and HCM frame as functions of M e g and with 
the measurement as a function of 2 ■ . The predictions from Ariadne are also 
shown. 
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Figure 11: Mean charged multiplicity, (n^) , in the current region of the Breit 
frame multiplied by 2 as a function of 2 ■ E% and in the current region of the HCM 
frame multiplied by 2 as a function of W . The results of e + e~ [34-38] and fixed- 
target DIS experiments [8-10] are shown. The factor 1.08 was estimated using MC 
predictions to correct the fixed-target data for the decay products of Kg and A. The 
predictions of Ariadne and Herwig are also shown. 
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